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Abstract — The use of miniaturized biomedical devices im- 
planted in the human body and wirelessly internetworked is 
promising a significant leap forward in medical treatment of 
many pervasive diseases. Recognizing the well-understood lim- 
itations of traditional radio-frequency wireless communications 
in interconnecting devices within the human body, in this paper 
we propose for the first time to develop network protocols for 
implantable devices based on ultrasonic transmissions. We start 
off by assessing the feasibility of using ultrasonic propagation in 
human body tissues and by deriving an accurate channel model 
for ultrasonic intra-body communications. Then, we propose a 
new ultrasonic transmission and multiple access technique, which 
we refer to as Ultrasonic WideBand (UsWB). UsWB is based 
on the idea of transmitting information bits spread over very 
short pulses following a time-hopping pattern. The short impulse 
duration results in limited reflection and scattering effects, and its 
low duty cycle reduces the thermal and mechanical effects, which 
are detrimental for human health. We then develop a multiple 
access technique with distributed control to enable efficient 
simultaneous access by interfering devices based on minimal 
and localized information exchange and on measurements at 
the receiver only. Finally, we demonstrate the performance of 
UsWB through a multi-scale simulator that models the proposed 
communication system at the acoustic wave level, at the physical 
(bit) level, and at the network (packet) level. 

I. Introduction 

Body area networks have received considerable attention 
in the last five years Q, Q, driven by the fascinating 
promise of a future where carefully engineered miniaturized 
biomedical devices implanted, ingested or worn by humans 
can be wirelessly internetworked to collect diagnostic infor- 
mation (e.g., measuring the level of glucose in the blood of 
diabetic patients) and to fine-tune medical treatments (e.g., 
adaptively regulate the dose of insulin administered) over 
extended periods of time l3l. Bl. 

Yet, most research to date has focused on communica- 
tions along the body surface among devices using traditional 
electromagnetic radio-frequency (RF) carrier waves, while 
the underlying challenge of enabling networked intra-body 
miniaturized (at the micro or nano scale) sensors and actu- 
ators that communicate through body tissues is substantially 
unaddressed. The main obstacle to enabling this vision of net- 
worked implantable devices is posed by the physical nature of 
propagation in the human body, which is composed primarily 
(65%) of water, a medium through which RF electromagnetic 
waves notoriously hardly propagate, even at relatively low fre- 
quencies. Accordingly, most research has focused on reducing 
the radiated power to avoid overheating of tissues 0. 

These formidable challenges cannot be overcome unless 
a major paradigm shift in networking through body tissues 
is made to address the limitations of RF propagation in the 
human body. Therefore, in this paper we present the first 



study proposing the use of ultrasonic waves to wirelessly 
internetwork in-body devices. 

Acoustic waves, typically generated through piezoelectric 
materials, are known to propagate better than their RF coun- 
terpart in media composed mainly of water. Since World 
War II, piezoelectrically generated acoustic waves have found 
application, among others, in underwater communications 
(typically at frequencies between and 100 kHz O), in 
indoor localization in sensor networks Q, and, massively, in 
ultrasonic medical imaging [8]. While communication at low 
frequencies requires sizable transducers, innovations in piezo- 
electric materials and fabrication methods, primarily driven 
by the need for resolution in medical imaging, have made 
miniaturized transducers, at the micro 0, |[T0l and even nano 
scales I'll] a reality. Moreover, the medical experience of the 
last decades has demonstrated that ultrasound is fundamentally 
safe, as long as acoustic power dissipation in tissue is limited 
to specific safety levels El, 0. It is also known that ultrasonic 
wave heat dissipation in tissues is minimal compared to RF 
waves fT2l . 

We lay our foundation on the fundamental physics of 
ultrasound propagation and move up layers of the protocol 
stack with the following core contributions: 

— Feasibility of Ultrasonic Communications in the Human 
Body. We start off by discussing the fundamentals of ultrasonic 
propagation in tissues, and by exploring important tradeoffs, 
including the choice of a transmission frequency, transmission 
power, bandwidth, and transducer size. 

— Ultrasonic Channel Modeling. We derive a model of 
the channel impulse respons for ultrasonic communications 
in tissues and observe that the inhomogeneity of the human 
body, characterized by a multitude of very small organs, causes 
severe multipath effect, thus making detection and decoding 
a challenging operation. In addition, we observe that the low 
speed of sound in tissues leads to high delays that need to be 
considered in system design. 

— Ultrasonic Wideband Design. To address the severe 
effect of multipath, we design and propose a new ultrasonic 
transmission and multiple access technique, which we refer 
to as Ultrasonic WideBand (UsWB). Ultrasonic wideband is 
based on the idea of transmitting very short pulses following 
an adaptive time-hopping pattern. The short impulse duration 
results in limited reflection and scattering effects, and its low 
duty cycle reduces thermal and mechanical effects that are 
detrimental for human health. A variable length spreading code 
is superimposed to the time-hopping pattern to further combat 
the effect of multipath and to introduce waveform diversity 
among interfering nodes. 

— Ultrasonic Wideband Adaptation and Multiple Ac- 
cess Control. We then study and develop multiple access 
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techniques with distributed control to enable multiple access 
among interfering implantable devices. Our proposed scheme 
is based on the idea of regulating the data rate of each 
transmitter to adapt to the current level of interference by 
distributively optimizing the code length and duration of the 
time hopping frame. 

— System Performance Evaluation. We evaluate the pro- 
posed scheme through a multi-scale simulator that evaluates 
our system at three different levels, i.e., (i) at the wave level by 
modeling ultrasonic propagation through reflectors and scat- 
terers, (ii) at the bit level by simulating in detail the proposed 
ultrasonic transmission scheme, (iii) at the packet level by 
simulating networked operations and distributed control and 
adaptation to evaluate metrics such as network throughput and 
packet drop rate. 

The rest of the paper is organized as follows. In Sections 
[IT] and III we discuss fundamental aspects of ultrasonic phys- 
ical propagation in human tissues and channel modeling. In 
Section IV we outline the design of UsWB and in Section 
|V| we illustrate the proposed medium access control protocol. 
Performance evaluation results are discussed in Section IvH 



Finally, in Section VII we conclude the paper. 



II. Ultrasonic Propagation in Human Tissues 

Ultrasonic waves originate from the propagation of mechan- 
ical vibrations of particles in an elastic medium at frequencies 
above the upper limit for human hearing, i.e., 20 kHz. Even if 
each particle oscillates around its rest position, the vibration 
energy propagates as a wave traveling from particle to particle 
through the material. Acoustic propagation through a medium 
is governed by the acoustic wave equation (referred to as the 
Helmhotz equation), which describes pressure variation over 
the three dimensions, V 2 P — ^ = 0, where P(x, y, z, t) 
represents the acoustic pressure scalar field in space and time, 
and c is the propagation speed in the medium. 

When ultrasonic waves propagate through an absorbing 
medium, the initial pressure, Pq reduces to P(d) at a distance 
d E3 as P(d) = P e~ ad , where a ([Np • cm" 1 ]) is the 
amplitude attenuation coefficient that captures energy dissipa- 
tion from the ultrasonic beam and is a function of the carrier 
frequency / as a = af b fl3j . where a ([Np • m _1 MHz _b ]) 
and b are tissue attenuation parameters. 

TABLE I: Frequency Limits for A=100 dB 



Communication Range 


Distance 


Frequency Limit 


Short Range 


/iin - mm 


> 1GHz 


Medium Range 


mm - cm 


~ 100MHz 


Long Range 


>cm 


~ 10MHz 



In we observed that attenuation can be significant and it 
increases with the distance between transmitter and receiver. 
Moreover the beam spread is inversely proportional to the 
ratio between the diameter of the radiating surface and the 
wavelength of the operating frequency 0. Consequently, 
since most biomedical sensing applications require directional 
transducers, one needs to operate at high frequencies to 
keep the transducer size small. However, as expected, higher 
frequencies lead to higher attenuation. 

In Table [I] we summarize our findings in on the 
maximum "allowed" carrier frequency for a 100 dB maximum 
tolerable attenuation. For distances ranging between some fim 
up to a few mm (i.e., short range communications) frequencies 
higher than 1 GHz are allowed. When distances are higher 



than 1 mm but still lower than some cm, i.e., for medium 
range communications, the transmission frequency should be 
decreased to approximately 100 MHz. For distances higher 
than a few cm, i.e., for long range communications, the 
transmission frequency should not exceed 10 MHz. Our results 
in are independently confirmed in 0. 

III. Channel Modeling 

To characterize the ultrasonic communication channel in tis- 
sues - for which, unfortunately, there is basically no literature 
available to date, except for [14] - we derive a deterministic 
channel model based on wave propagation theory. Propagation 
of acoustic waves through biological tissues is governed 
by three coupled first-order equations, i.e., the continuity 
equation, the force equation and the equation of state fT5lL 
which represent relationships among acoustic pressure, P, 
acoustic particle velocity u, and medium density p, and can 
be rearranged to obtain the Helmhotz equation in Section [TTJ 
The solution represents the acoustic field in time, evaluated at 
each spatial coordinate of the propagation medium. 

Traditionally, partial differential equations are solved using 
numerical methods such the finite-difference-method (FDM). 
We take a different, computationally more efficient approach 
based on the pseudo- spectral and k-space methods 1 16]. Basi- 
cally, the pseudo-spectral (PS) method reduces computational 
complexity in the spatial domain by using Fourier series 
expansions and FFTs, while the k-space method operates 
in the time domain by using k-space propagator functions 
(instead of classical finite differences) to approximate temporal 
derivatives. As a consequence, larger time steps can be used, 
reducing the simulation time with controllable accuracy. A 
powerful tool implementing this method is k-Wave ifTTl . 

We modeled a reference ultrasonic propagation channel in 
the human body to derive an accurate characterization of 
the channel impulse response, which was then used in PHY 
and MAC layer simulation studies. Specifically, we modeled 
a section of the human arm, including bones, muscles, fat 
and skin. We considered an heterogeneous two-dimensional 
rectangular area of length 20 cm and width 10 cm, where 
dens ity and sound velocity are distributed as shown in Fig. 
| l(a)[ The bone half-section has a 18 mm width, the muscle 
half-section has a 22mm width, while fat and skin have a half- 
section of 7 mm and 3 mm width, respectively. The medium 
parameters, sound speed c, density p and attenuation a and b 
are as in fT8ll . 

We obtained the channel impulse response by releasing 
an ideal Dirac pulse in the left top corner of the muscle 
section. The receiver is located in the left bottom corner of 
the muscle section. Transmitter and receiver are located 20 cm 



away from one another. In Fig. 1(b) we show snapshots of the 
acoustic field time propagation in the modeled environment. 
We observe that the effect of the bone and tissues is to partially 
reflect and scatter the acoustic wave transmitted by the source. 
We then obtain the channel impulse response by recording the 
time series of the signal at the receiver sensor. 

Figure 1(c) reports the resulting impulse response in this 
scenario. We observe that the effect of multipath and scattering 
is to introduce attenuated signal replicas spaced in time. 
Because of the very short duration of the transmitted pulses, 
replicas do not interact destructively, which provides a strong 
motivation for the proposed transmission scheme presented 
in Section [TV] We can model the channel response as a 
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Fig. 1: (a) Sound speed and density distributions in a simulated arm. (b) Snapshots of the time evolution of acoustic field, (c) Resulting ultrasonic 
channel impulse response in the arm. 



complex- valued low-pass equivalent impulse response. We can 
then characterize the channel through the mean excess delay 
(r m ) and the RMS delay spread (trms)- For the channel 
simulated above, we obtained r m = 1.2779 • 10 -5 s and 
trms — 2.6883 • 10 _5 s. Since the coherence bandwidth of 
the channel is proportional to the inverse of trms, we should 
consider the above channel as frequency selective for signals 
of bandwidth above approximately 7 kHz. 

IV. Ultrasonic Wideband 

Based on the considerations summarized in Sections [n] and 
III we design and propose a new ultrasonic transmission and 
multiple access technique, which we refer to as Ultrasonic 
WideBand (UsWB). The key design objectives of UsWB are 
(i) to enable low-complexity and reliable communications in 
ultrasonic channels against the effect of multipath reflections 
within the human body; (ii) to limit the thermal effect of 
communications, which is detrimental to human health; (iii) to 
enable distributed medium access control and rate adaptation 
to combat the effect of interference from co-located and 
simultaneously transmitting devices. UsWB is jointly designed 
to provide physical layer functionalities and medium access 
control arbitration and adaptation to enable multiple concur- 
rent co-located transmissions with minimal coordination. 

Ultrasonic Pulsed Transmissions. Ultrasonic wideband is 
based on the idea of transmitting very short ultrasonic pulses 
following an adaptive time-hopping pattern together with a 
superimposed adaptive spreading code. Baseband pulsed trans- 
missions enable high data rate, low-power communications, 
low-cost transceivers, and have been proposed for RF short- 
range, high data rate communications (T9J, [20], although 
with much shorter pulse durations (and consequently larger 
bandwidth) than achievable in ultrasonic communications. 

The characteristics of pulsed transmissions appear to ideally 
address the requirements discussed above. Their fine delay- 
resolution properties are well- suited for propagation in the 
human body, where inhomogeneity in terms of density and 
propagation speed, as well as the pervasive presence of very 
small organs and particles, cause dense multipath and scatter- 
ing. When replicas of pulses reflected or scattered are received 
with a differential delay at least equal to the pulse width, they 
do not overlap in time with the original pulse. Therefore, for 
pulse durations in the order of hundreds of nanoseconds 0, 
pulse overlaps in time are reduced and multiple propagation 
paths can be efficiently resolved and combined at the receiver 
to reduce the bit error rate. Also, the low duty cycle of pulsed 



transmissions reduces the impact of thermal and mechanical 
effects, which are detrimental for human health |2H . Finally, 
carefully designed interference mitigation techniques may 
enable MAC protocols that do not require mutual temporal 
exclusion between different transmitters. This is crucial in the 
ultrasonic transmission medium since the propagation delay 
is five orders of magnitude higher than in RF in-air channels 
and carrier- sense-based medium access control protocols are 
ineffective 1 6 ] . In addition, data rate can be flexibly traded for 
power spectral density and multipath performance. 

A. Physical Layer Model 

Adaptive Time-Hopping. Consider a slotted time divided 
in chips of duration T c , with chips organized in frames of 
duration Tf = Nh • T c , where Nh is the number of chips 
per frame. Each user transmits one pulse in one chip per 
frame, and determines in which chip to transmit based on a 
pseudo-random time hopping sequence (THS), i.e., a sequence 
generated by seeding a random number generator with the 
user's unique ID. 

The train of pulses is modulated based on pulse position 
modulation (PPM), i.e., a '1' symbol is carried by a pulse 
delayed by a time S with respect to the beginning of the chip, 
while a '0' symbol begins with the chip. The signal s( k \t,i) 
generated by the k th user to convey the i th symbol is expressed 
as 

= p{t - cf } T C - iT f - dfh), (1) 



where p(t) is the second derivative of a Gaussian Pulse, {c-^} 
is the time hopping sequence of the k th source, with < 
c-^ < Nh — 1, and {d[ k ^} is the information-bearing sequence, 
d^ G {0, 1}. The resulting data rate, in pulses per second, 
is expressed as R(Nh) = = N ^ T • By regulating the TH 
frame length Nh, i.e., the average inter-pulse time, a user can 
adapt its transmission rate, and as a consequence modify the 
average radiated power and therefore the level of interference 
generated to other ongoing communications. We observe at 
this point that an individual user has little incentive to increase 
its frame size, since that results in a lower achievable data rate, 
without any major benefit for the user itself (since the level of 
interference perceived depends primarily on the frame length 
of the other users, and not on its own). However, a longer time 
frame reduces the interference generated to the other users. 
Therefore, selfish/greedy frame adaptation strategies do not 
work well in this context - cooperative strategies are needed. 
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Adaptive Channel Coding. An adaptive channel code 
l22ll can further reduce the effect of mutual interference 
from co-located devices, by dynamically regulating the coding 
rate to adapt to channel conditions and interference level. 
Various channel coding solutions have been proposed fT9lL 
l22l with different performance levels and computational 
complexity. We rely on simple pseudo- orthogonal spreading 
codes because of their excellent multiple access performance, 
limited computational complexity, and inherent resilience to 
multipath. Each symbol (i.e., bit) is spread by multiplying it 
by a pseudorandom code before transmission. At the receiver 
side, with prior knowledge of the code used at the transmitter, 
the signal can be de- spread, and the original information 
recovered. With different orthogonal (or pseudo-orthogonal) 
codes, multiple nodes can transmit simultaneously on the 
same portion of the spectrum, with reduced interference. We 
explored two alternative modulation schemes, which we refer 
to as PPM-BPSK-spread and PPM-PPM- spread. 

In PPM-BPSK-spread, the binary spreading code at the k th 
node, {a^}, is defined as a pseudorandom code of N s chips 
with cij G { — 1,1}. Accordingly, the information bit is spread 
using BPSK modulated chips, and by combining with time 
hopping, ([T]) can be rewritten as 

N a -1 

*<*>(M)= E afp{t-cfT c -jT f -dfh), (2) 
3=0 

where 5 is the PPM displacement of a pulse representing a ' 1 ' 
bit, while chip information is carried in the pulse polarity. 

In PPM-PPM-spread, the information bit is spread using 
PPM-modulated chips. In this case, the binary spreading code 
can be defined as a pseudorandom code of N s chips with 
(ij G {0, 1}. With time hopping, |T|) can be rewritten as 



(M) 



N s -1 
3=0 



C 3 ±c 



O*) 0) 



where © represents a XOR operation. 

In Fig. [2] we show an example of a combined time hopping 
and PPM-BPSK-spread coding strategy. Since the spreading 
operation associates N s chips to one information bit, the 
information rate will be further reduced by a factor N s , i.e., 



R(N h ,N s ) 



1 



1 



N a T f N s N h T c ' 



(4) 



while the energy per bit is increased by a factor N s . Note 
that there is a tradeoff between robustness to multi-user 
interference (which increases with longer spreading codes), 
and energy consumption and information rate. In Section [V] 
we discuss joint dynamic adaptation of frame and code length. 
With BPSK-modulated chips the absolute received phase 
information is needed for decoding. Therefore, a coherent 
receiver with accurate channel knowledge is needed. Instead, 
with a pure PPM-modulated signal, a simple non-coherent 
energy detector receiver is sufficient. . The latter requires 
frame synchronization only, and its hardware complexity is 
significantly lower 1 231 . 

Signal to Interference-plus-Noise Ratio. We can express 
the signal to interference-plus-noise ratio (SINR) for impulsive 
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h = Tc 

Fig. 2: Example of two ongoing transmissions using time hopping 
sequences TH X = {3,12,8,2} and TH 2 = {11,3,12,7} and PPM- 
BPSK spreading codes Sd = {1, -1, -1, 1} and SC 2 = {1, 1, 1, -1}. 



transmissions at the receiver of link i as fT9l 
SINRi(N h , N s ) = N s PiSid NhaTc 



(5) 



where Pi is the average power per pulse period emitted by the 
i th transmitter, gij is the path gain between the i th transmitter 
to the j th receiver, 77 represents background noise energy and 
a 2 is an adimensional parameter depending on the shape of 
the transmitted pulse and the receiver correlator. The set Z$ 
represents the set of links whose transmitter interferers with 
the receiver of link i. Remind that Tf — NhT c is the frame 
duration. The expression in §5§ depends on the vectors of frame 
and code lengths of all the ongoing communication in the 



network, i.e., Nh,N s . The term 



accounts for the level 



of interference generated by each interferer k to the receiver 
of link i, i.e., the number of pulses transmitted by the k th 
transmitter during the time frame of the i th user. Note that 
we will not consider power control strategies in our treatment, 
since 1 20 ] showed that in the linear regime, when the objective 
is to maximize the aggregate data rate, the optimal solution 
always corresponds to points where individual devices transmit 
at the maximum power, or do not transmit at all. 

Note that increasing (or decreasing) the spreading code of 
the node of interest, N s ^, leads to an increase (decrease) in 
the SINR. When the node of interest increases (decreases) its 
frame length, Nh,i, while the other nodes do not, we expect 
no variation in the SINR (there is in fact a slight increase 
(decrease) in the SINR, which can be neglected under high 
SNR conditions, r\ <C Ylkex- PkQki)- Finally, when the frame 
length of the interfering nodes is increased (decreased), the 
SINR increases (decreases). 

V. MAC and Rate Adaptation 

In this section, we discuss UsWB medium access control 
principles and rate adaptation. Based on our discussion so 
far, there is a tradeoff between (i) resilience to interference 
and channel errors, (ii) achievable information rate, and (iii) 
energy efficiency. In this section, we introduce medium access 
control and rate adaptation strategies designed to find optimal 
operating points along efficiency-reliability tradeoffs. We first 



consider rate-maximizing adaptation strategies in Section V-A 



Then, we propo se two different energy-minimizing strategies 
in Section IV-BI 

A. Distributed Rate -Maximizing Adaptation 

The objective of the rate- adaptation algorithm under con- 
sideration is to let each active communication maximize its 
transmission rate by selecting a pair of code and frame lengths, 
based on the current level of interference and channel quality 
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measured at the receiver and on the level of interference to the 
other ongoing communications generated by the transmitter. 
We consider a decentralized ultrasonic intra-body area net- 
work. Denote by M the set of \M\ existing connections, and 
by Nh,max and N Sy7nax the maximum lengths supported, 

< N hti < N htmax , Vi e A/\ N h e N, (6) 

< N 3j i < N s , max , \/i G A/", N s g N, (7) 

where N is the set of natural numbers. According to the 
transmission scheme discussed in Section |IV-A| each node 
transmits at a rate R as in ^ and each receiver experiences 
an SINR as in ([5]). Each node has a minimum data rate 
requirement, i.e., Ri(N^ ji? N 3j i) > R m in, and a minimum 
SINR requirement, i.e., 5JiV^(N h ,i, N s ,i) > SINR min . 

Explicitly Cooperative Problem. Denote the frame length 
and the code length selected by the receiver of the communi- 
cation r as Nh, r and N s ^ r . 

The objective of each user is to locally optimize the infor- 
mation rate of the communication by solving the following 
problem: 



maximize R r (Nh r , N s r ) 

N h , r ,N s , r ' 



(8) 



subject to R r (N h , r ,N s , r ) > R min (9) 
SINR r (N h , r , N Syr ) > SINR min (10) 
SINRi(N hir ,N ai r) > SINR min Vi e X r (11) 

where X r is the set of the connections interfering with the r th 
connection. The constraints on the maximum frame and code 
length in ^ and ^ are also implicitly considered. We refer 
to this as the explicitly cooperative problem. 

Implicitly Cooperative Problem. If all nodes measure the 
same level of interference, that is, all network nodes are 
close enough to be all in the same transmission range, and 
all nodes have the same minimum rate and minimum SINR 
requirements, the level of interference that can be tolerated by 
each receiver is the same. Therefore, this information does not 
need to be exchanged. The problem in ([5]) becomes 

maximize R r (Nh r , N s r ) 

N h , r ,N s , r ' 

subject to R r (N h , r , N s , r ) > R min (12) 
SINR r (N Kr ,N s , r ) > SINR min (13) 

The system of SINR inequality constraints in ^ becomes 
here a single inequality constraint. The new problem can be 
interpreted as finding the optimal pair of code and frame length 
that maximize the rate, given a minimum rate and a minimum 
SINR, under the assumption that all the other nodes will be 
acting in the same way. As we will show in Section [Vl| since 



the problem to be solved is the same for each node, a globally 
optimal pair of code and frame lengths will be found. We refer 
to this as the implicitly cooperative problem. 

We now provide a detailed analytical description of the 
explicitly cooperative problem. The implicitly cooperative 
problem can then be derived as a special case. 

The explicitly cooperative problem as stated in ([5]) is an 
integer program, i.e., variables Nh and N s are integer. If 
the domain of the problem is small, it can be solved by 
enumeration, i.e., trying all the possible combinations of Nh 
and N s . When the domain of the problem increases in size 
a relaxation method can be used to transform the integer 



problem into a fractional problem. The relaxation operation 
consists of replacing the constraints in ^ and ^ with 

< N h < N h 

,max i < N s < N s 

,max i 

N s ,N h eR, (14) 

where R is the set of real numbers. Note that the optimal 
solution to the relaxed problem is not necessarily integral. 
However, since the feasible set of the relaxed problem is 
larger than the feasible set of the original integer program, 
the optimal value of the former, p* lx is a lower bound on the 
optimal value of the latter, p\ nt , i.e., L = p* lx < p* nt . The 
relaxed solution can be used to guess the integer solution by 
rounding its entries based on a threshold t G [0,1]. 

If the rounded solution, x* nt , is feasible for the original 
problem, i.e., all the constraints are satisfied, then it can be 
considered a guess at a good, if not optimal, point for the 
original problem. Moreover, the objective function evaluated 
at p* nt is an upper bound on p* nt , U = p* nt > p* nt Note that 
x \ nt is cannot be more than (U-L)-suboptimal for the original 
problem. We can express the relaxed problem as a geometric 
program (24), that is, minimizing a posynomial function under 
posynomial inequality constraints. First, the problem must be 
restated in terms of minimizing the inverse of the data rate, that 
is, R~ 1 (N} ly r^ N Syr ). Then, the new objective function and all 
the constraints can be expressed in monomial and posynomial 
form. The objective function, as well as the minimum rate 
function, is a monomial function, N^l N~* . The same holds 
for the constraint functions defined is ( fT4] ), and for the SINR 
functions defined in ([5]). 

Let us first consider the SINR constraint in ( fTO] ), and note 
that in the distributed optimization problem the only variables 
are the code and the frame length of the receiver of interest, 
N s ^ r and Nh, r - If we define 

a r = P r g r ^ r T Cj f3 k = cr 2 T c P k g ryk /N hykj 

then {TO} can be rewritten as a posynomial constraint 

VN-^N^I + AT"* J2 & < • (15) 

kex r 



SINR m ir 



For the remaining (\X r \ — 1) SINR constraints, if we define 

7* Pi9i,iT c N hii N Sii ei cr 2 T c P r gi^N h)i 

Si = E k ez l a2 TcPkgi,kNh,i/Nh y k (16) 
then (T\) can be rewritten as 

•q + Si + eiN^l _ 1 



< 



(17) 



7i SINR min 

which is a linear, and thus monomial, constraint. In particular, 
this constraint shows that the frame length variation only leads 
to an increase (decrease) in the level of interference produced. 
Moreover, the constraint suggests that the optimal frame length 
is upper bounded by the co-located node that is experiencing 
the highest level of interference. Based on these observations, 
the cooperative optimization problem becomes 

(18) 



minimize R r 1 (Nh r , N s r ) 

N h , r ,N s , r r 

subject to R-\N h , r ,N s , r ) < R m \ n 

vN-}N-l+N-^^< 



a r 



kex r 



SINRn 



N h , r > 



7z SINR 



i 

min 



TJ-Si 



(19) 
(20) 

(21) 
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Finally, a geometric program can be in general transformed 
into a convex program through a logarithmic transformation of 
the optimization variables, which can_be solved in polynomial 
time through interior-point methods 



B. Distributed Energy -minimizing Rate Adaptation 

We now concentrate on rate adaptation with the objective 
of reducing the energy consumption of UsWB. We define (i) 
E p , the energy per pulse; (ii) E^, the energy per bit, i.e., = 
E p • N s ; and (iii) E s , the average energy emitted per second, 
i.e., E s =E v /(T c .N h ). 

The energy per bit is a linear function of the spreading code 
length, therefore of the number of pulses transmitted per each 
bit. The average power emitted per second is a function of 
the inverse of the frame length, hence of the number of pulses 
transmitted per second. However, both depend on the value of 
Ep, the energy transmitted per pulse, which is related to the 
(electrical) power absorption of the ultrasonic transducer. 

Modeling the Energy Consumption of a Piezoelectric 
Transducer. We consider electronically-driven piezoelectric 
transducers. Since their electronics are known to have minimal 
current leakage, the majority of power consumption comes 
from the piezoelectric element. Outside the region of reso- 
nance, a piezoelectric ceramic transducer can be viewed (from 
the electrical point of view) as a parallel plate capacitor with 
capacitance Cq. Thus, the main source of power consumption 
comes from charging such capacitor (26). Then, ignoring 
charge and discharge losses, and considering a capacitor with 
voltage supply V and pulse repetition frequency / (which 
corresponds to the charge and discharge frequency of the 
capacitor), the power consumption P c can be expressed as 



P c = fC V 2 . 



(22) 



The static capacitance and the voltage supply values should 
be based on transducer- specific considerations. The static 
capacitance value of a disc- shaped transducer is given as 
Co = Ae t ° K [F], where eo is the permittivity in free air 
(8.8542 F/W), K is the dielectric constant of the material 
(adimensional), A [m 2 ] is the area of the disc, and th [m] is 
its thickness. The dielectric constant K depends on material, 
frequency and mechanical state of the transducer. 

We derive appropriate limits for the voltage supply based on 
safety concerns, which impose a limit on the radiated acoustic 
power. As reported in El. no tissue damage occurs in intra- 
body ultrasonic propagation as long as the acoustic power 
dissipation in tissue is limited to 10 4 W/m 2 . From this limit, 
we derive the corresponding maximum pressure magnitude 
that can be radiated by the transducer, and consequently 
the maximum voltage input. For example, if we consider 
arm muscle as the primary propagation medium, through 
the quadratic relation between the acoustic intensity and the 
acoustic pressure, I = (Prms) 2 /pc, along with density and 
speed of sound parameters reported in (TBI , the maximum 
pressure magnitude that can be radiated by the transducer is 
found to be approximately 0.13 MPa. 

We can derive the related transducer voltage input, cor- 
responding to the maximum radiated pressure, through the 
constitutive equation of piezoelectric materials. The latter 
expresses the relationship among mechanical strain and elec- 
trical displacement for the piezoelectric element considering 
the electrical and mechanical stress and is usually expressed 



in tensor notation l27l . Here, we are interested in the so- 
called converse effect, i.e., the material strain proportional 
to an applied voltage. If we consider a symmetric crystal 
structure, the constitutive relations reduce to a few parameters. 
In particular, the so-called #33 parameter [Vm/N] represents 
the piezoelectric voltage coefficient, when the polarization 
field and the piezoelectrically induced strain are both parallel 
to the disc axis (usually referred to as the third axis). Based 
on this, one can express the electric field along the third axis 
E [V/m] as E = g 33 P, where P again represents the output 
pressure of the transducer. Since the electric field E can in 
turn be expressed as the ratio between voltage and electrode 
distance, we have V = g^Pth- Finally, from (22), the energy 
consumed to generate a single pulse is obtained as 



E P = C (gssPh 



[J]. 



(23) 



Energy-minimizing Rate Adaptation. Based on this 
model, we design a rate adaptation strategy where the objective 
of r is to minimize (i) the energy per bit, E^, or (ii) the average 
energy emitted per second E s . The problem consists of finding 
the optimal frame length and the optimal spreading code 
length N s that minimize (E s ) while meeting the minimum 
SINR constraints and keeping the data rate over a given 
threshold. The energy problem is formally expressed below. 

minimize Eb(N s r ) (or E s (Nh r )) 

N h ,r,N s , r 

subject to ([9]> ([TOj) <[TT|) 

Similarly, the problem above can be relaxed to a geometric 
program as discussed in Section |V-A 



C. Medium Access Control Protocol 

In UsWB, distributed medium access control coordination 
is achieved by exchanging information on logical control 
channels, while data packets are transmitted over logical data 
channels. We consider unicast transmissions between a trans- 
mitter Tx and a receiver Rx. When Tx needs to transmit a 
packet, it first needs to reserve a dedicated channel to Rx. The 
connection is opened through the common control channel, 
which is implemented through a unique TH-sequence and a 
spreading code known and shared by all network devices. 

In the two-way handshake procedure, T x sends a 
Request-to-Transmit (R2T) packet to R x , which 
contains its own ID. If R x is idle, a Clear-to-Transmit 
( C 2 T ) control packet is sent back to T x . In case of failure and, 
thus, timer expiration, T x will attempt a new transmission after 
a random backoff time, for a maximum of Nr times. During 
these initial steps, since an estimate of the current interference 
is not available, the transmitter transmits at the minimum 
data rate, conservatively using a maximum frame length and 
spreading code length, thus generating low interference. After 
receiving the C2T, the transmitter switches to a dedicated 
channel by computing its own TH sequence and spreading 
code obtained by seeding a pseudo random sequence generator 
with its ID. As a consequence, both T x and R x leave the 
common channel and switch to a dedicated channel. Once 
the connection has been established, T x sends the first packet 
using maximum frame and spreading code length. The receiver 
R x computes the opti mal fr ame and spreading code lengths 
as discussed in Section [TV- A| This information is piggybacked 
into ACK or NACK packets. 
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In the explicitly cooperative case, discussed in Section |V-A 
once the communication has been established, R x does not 
leave the common control channel. Instead, it keeps "listening" 
to both the dedicated and common control channel at the 
same time. In the dedicated control channel, R x sends to T x 
the optimal frame and code lengths to be used for the next 
transmission. In the common control channel, R x exchanges 
with other co-located receivers information on the level of 
tolerable interference. 

VI. Performance Evaluation 

In this section we evaluate the proposed system through a 
custom-designed multi- scale simulator that considers UsWB 
performance at three different levels, i.e., (i) at the acoustic 
wave level by modeling ultrasonic propagation in tissues 
through reflectors and scatterers, (ii) at the bit level by 
simulating in detail the ultrasonic transmission scheme, (iii) 
at the packet level by simulating networked operations and 
distributed medium access control and adaptation. 

The acoustic wave level simulation is performed as de- 
scribed in Section [III] A channel impulse response is obtained 
by simulating propagation in the human arm. Transmission 
at the bit level is modeled through a custom physical layer 
simulator, which allows to obtain an empirical model of 
the Bit Error Rate (BER) against different values of the 
TH frame length and spreading code length, for different 
levels of interference. The physical layer simulation models a 
transmitter and a receiver (located 20 cm apart) communicating 
over an acoustic channel through UsWB. Each sender node 
transmits information as discussed in Section HV] Simulations 
are performed to obtain an estimate of the achievable BER 
upon varying the frame length and the spreading code length 
with a different number of interferers transmitting on the same 
chann el. We implement both solutions presented in Section 
IV-A[ i.e., PPM-B PS K- spread with coherent receiver and the 



PPM-PPM- spread with non-coherent receiver. An extensive 
simulation campaign was conducted and BER values were 
obtained as a function of the spreading code and frame lengths 
for different number of simultaneously active connections. 
Figure [3] reports the BER values when 4 pairs of nodes are 
transmitting at the same time, and as expected the coherent 
PPM-BPSK-spread outperforms the non-coherent PPM-PPM- 
spread. The empirical model of the BER as a function of 
frame length and spreading code length for a given level 
of interference is then imported in a Java-based event-driven 
packet-level simulator, which models all aspects of the UsWB 
communication protocol stack. 

Network Simulation Topology. We considered two dif- 
ferent settings for the network level simulations. First, we 
consider a 2-D topology with 18 static nodes randomly lo- 
cated inside a square of side 20 cm. A maximum of nine 
communicating pairs are considered. We also assume that 
the transmission range is greater than the maximum distance 
between the nodes, thus all nodes measure the same number 
of interferers. The second setting consists of three 2-D squared 
clusters of nodes with side 10 cm, displaced 20 cm apart 
from each other (distance center- to-center). In each cluster, 
the nodes are randomly deployed according to a gaussian 
distribution. The intermediate located cluster (#1) contains one 
communicating pair, while the other clusters, (#2) and (#3), 
contain a variable number of communicating pairs. We also 
assume the transmission range to be of 30 cm, thus all the 
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Fig. 3: (a) BER vs. TH frame length and (b) BER vs. spreading code 
length for PPM-BPSK spread and PPM-PPM spread. 

nodes of adjacent clusters interfere with each other, while 
nodes of non adjacent clusters do not. Under this assumption, 
communicating pairs from different clusters measure different 
levels of interference. 

During the simulations, the frame length and spreading 
code length are adapted to the interference level measured 
at the rece iver a ccording to one of the objectives discussed 
in Section |V-A| In particular, the first setting matches the 
condition required by the implicitly cooperative problem in 
( [T6] ), while for the second setting the explicitly cooperative 
problem in ^ is considered. The simulation time is set to 
100 s and nodes start establishing connections at a random 
time instant but no later than 2 s after the simulation start time. 
We consider an infinite arrival rate at each transmitter, i.e., 
transmitters are always backlogged. The maximum allowed 
frame and spreading code length is set to 15 slots and 20 
chips, respectively. The maximum supported rate, achieved 
when frame and spreading code length are both set to one, is 
equal to 2Mbit/s. The minimum SINR constraint is equivalent 
to a maximum BER constraint set to 10 -6 . 

A. Rate Maximization 

Implicitly cooperative solution. In Fig. [4] we plot and 
compare network throughput and packet drop rate for the rate- 
maximizing strategy, when frame and code length are adap- 
tively reg ulated based on the implicitly cooperative problem 
in Section V-A The rate-maximizing solution is presented for 
both the transmission schemes discussed in Section |IV-A[ the 
PPM-BPSK-spread with coherent receiver and the PPM-PPM- 
spread with non-coherent receiver. As expected the coherent 
PPM-BPSK solution performs better in terms of throughput. 
This happens because the BER constraint is satisfied for lower 
values of frame length and code length, which leads to higher 
data rates according to ([?]). 

In terms of packet drop rate, for the coherent PPM-BPSK 
system the BER constraint is satisfied for any number of active 
connections considered. Instead, for the non-coherent PPM- 
PPM system, when the number of active connections is greater 
than seven the BER constraint cannot be satisfied anymore, 
and therefore the packet drop rate increases. Note that this 
problem can be overcome by simply relaxing the constraint 
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Fig. 5: Time evolution of TH frame and spreading code lengths for 
implicitly cooperative problem. 



on the maximum size of code length or frame length, without 
relaxing the constraint on the maximum BER. Clearly this 
leads to a lower data rate. 

We then focus on the analysis of the dynamic adaptation 
of the frame and code length performed distributively by each 
node. Focusing on the coherent system, in Fig. [5] we show the 
evolution in time of frame and code length when 9 different 
connections are asynchronously activated with a deterministic 
5 s delay between each other. As expected, each connection 
starts with the maximum supported frame and code length 
and then adaptively reaches the optimal value based on the 
interference level measured at the receiver. Since the number 
of interferers is the same for each receiver, the locally optimal 
solution is also globally optimal. 

Explicitly cooperative solution In Fig. [6] we plot the 
network throughput obtained when frame and code length 
are adaptively regulated based on the explicitly cooperative 
problem in ([5]). The solution is presented for the PPM-BPSK- 
spread transmission strategies discussed in Section |IV-A 



Throughput is evaluated by varying the number of active 
connections in four consecutive time steps. In particular, we 
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explicitly cooperative problem. 



assume that there is always one active connection in cluster 
#2. Cluster #1 and #3, activate a new connection at each 
time step. In Fig. [6] we also report the average number 
of interferers measured by the receivers in each cluster. As 
expected, the cluster located in the middle, hence located 
within the transmission range of both the other two clusters, 
measures a higher level of interference, and therefore achieves 
a lower throughput. However, since the BER constraint is 
satisfied by all nodes in the three clusters, a low packet drop 
rate is achieved. Finally, in Fig. [7] the dynamic behavior of 
the frame and code length adaptation is shown. We consider 
two new connections activated asynchronously every 5 s in 
cluster #1 and #3. In each cluster, the frame and code length 
are adapted according to the increasing level of interference 
measured in the channel. In particular, we observe the effect 
of the constraint in ( [2T] ), which forces the frame length to 
be greater or equal than the frame length of the connection 
that is experiencing the highest level of interferences, i.e., the 
connection in cluster #2. 
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Fig. 8: Throughput and Eb vs. number of active connections for rate 
optimized and ^-optimized solutions. 

B. Rate-optimal Vs. Energy-optimal Results 

In Figs. [8] we plot the throughput obtained by adapting the 
transmission rate acc ording to the energy-minimizing strategy 
introduced in Section V-B The energy-minimizing strategy is 
also co mpare d to the rate-maximizing strategy introduced in 
Section V-A in terms of data rate and energy consumption. 
The throughput achievable in case of the rate-maximizing 
solution is comparable to what is obtained with the energy- 
optimized solution. When the number of active connections, 
N, is lower than 2, the throughput of the ^-optimal strategy 
is close to the throughput of the rate-maximizing scheme. 
For N higher than 2, the two considered strategies show a 
similar throughput performance. Since the maximum BER 
constraints are always satisfied, the rate-maximizing and the 
energy-minimizing strategies lead to packet drop rates close to 
zero. Finally, the energy consumption obtained with the E^- 
optimal solutions, is always better that what can be obtained 
using the rate-maximizing scheme. 

VII. Conclusions 

We proposed a paradigm shift in networking through body 
tissues to address the limitations of RF propagation in the 
human body. We presented the first attempt at enabling 
networked intra-body communications among miniaturized 
sensors and actuators using ultrasonic waves. 

The contribution of our work is manyfold. We assessed the 
feasibility of using ultrasonic communications within the hu- 
man body; we derived an accurate channel model for ultrasonic 
communications in the human body and built on it to propose 
a new ultrasonic transmission and multiple access technique, 
denoted as UsWB, based on transmission of short duration 
pulses following a time hopping pattern. Simulation results 
demonstrate the high throughput, low packet drop rate, and low 
energy consumption of UsWB in intra-body communications. 
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